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NATIONAL  FOREWORD 


This  Indian  Standard  (Part  1/Sec  48)  which  is  identical  with  lEC  60793-1  -48  :  2007  'Optical  fibres  — 
Part  1-48:  Measurement  methods  and  test  procedures  —  Polarization  mode  dispersion'  issued  by 
the  International  Electrotechnical  Commission  (lEC)  was  adopted  by  the  Bureau  of  Indian  Standards 
on  the  recommendation  of  the  Fibre  Optics,  Fibres,  Cables  and  Devices  Sectional  Committee  and 
approval  of  the  Electronics  and  Information  Technology  Division  Council. 

The  text  of  lEC  Standard  has  been  approved  as  suitable  for  publication  as  an  Indian  Standard  without 
deviations.  Certain  conventions  are,  however,  not  identical  to  those  used  in  Indian  Standards. 
Attention  is  particularly  drawn  to  the  following: 

a)  Wherever  the  words  'International  Standard'  appear  referring  to  this  standard,  they  should  be 
read  as  'Indian  Standard'. 

b)  Comma  (,)  has  been  used  as  a  decimal  marker  while  in  Indian  Standards,  the  current  practice 
is  to  use  a  point  (.)  as  the  decimal  marker. 


In  this  adopted  standard,  reference  appears  to  certain  International  Standards  for  which  Indian 
Standards  also  exist.  The  corresponding  Indian  Standards  which  are  to  be  substituted  in  their  respective 
places  are  listed  below  along  with  their  degree  of  equivalence  for  the  editions  indicated: 


International  Standard 

lEC  60793-1-1  :2008  Optical  fibres 

—  Part  1  -1 :  Measurement  methods 
and  test  procedure  —  General  and 
guidance 

lEC  60793-1  -44  : 201 1  Optical  fibres 

—  Part  1  -44:  Measurement  methods 
and  test  procedures  —  Cut-off 
wavelength 


Corresponding  Indian  Standard 

IS/IEC  60793-1-1  :  2008  Optical 
fibres:  Part  1  Measurement  methods 
and  test  procedures.  Section  1 
General  and  guidance 

IS/IEC  60793-1-44  :  2011  Optical 
fibres  :  Part  1  Measurement  methods 
and  test  procedures.  Section  44  Cut- 
off wavelength 


Degree  of  Equivalence 
Identical 


do 


The  technical  committee  has  reviewed  the  provisions  of  the  following  International  Standards/Other 
Publications  referred  in  this  adopted  standard  and  has  decided  that  they  are  acceptable  for  use  in 
conjunction  with  this  standard: 


International  Standard/ 
Other  Publication 

lEC  60793-2-50 


lEC  60794-3 
lEC  61280-4-4 

lEC/TR  61 282-3 
lEC/TR  61 282-9 
lEC  61290-1 1-1 


Title 


Optical  fibres  —  Part  2-50:  Product  specifications  —  Sectional 
specification  for  class  B  single-mode  fibres 

Optical  fibre  cables  —  Part  3:  Sectional  specification  —  Outdoor 
cables 

Fibre  optic  communication  subsystem  test  procedures  —  Part  4-4: 
Cable  plants  and  links  —  Polarization  mode  dispersion 
measurement  for  installed  links 

Fibre  optic  communication  system  design  guides  —  Part  3: 
Calculation  of  link  polarization  mode  dispersion 

Fibre  optic  communication  system  design  guides  —  Part  9: 
Guidance  on  polarization  mode  dispersion  measurements  and  theory 

Optical  amplifier  test  methods  —  Part  11-1:  Polarization  mode 
dispersion  — Jones  matrix  eigenanalysis  method  (JME) 


(Continued  on  third  cover) 
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INTRODUCTION 

Polarization  mode  dispersion  (PIVID)  causes  an  optical  pulse  to  spread  in  the  time  domain. 
This  dispersion  could  impair  the  performance  of  a  telecommunications  system.  The  effect  can 
be  related  to  differential  phase  and  group  velocities  and  corresponding  arrival  times  5t  of 
different  polarization  components  of  the  signal.  For  a  sufficiently  narrow  band  source,  the 
effect  can  be  related  to  a  differential  group  delay  (DGD),  Ar,  between  pairs  of  orthogonally 
polarized  principal  states  of  polarization  (PSP)  at  a  given  wavelength.  For  broadband 
transmission,  the  delays  bifurcate  and  result  in  an  output  pulse  that  is  spread  out  in  the  time 
domain.  In  this  case,  the  spreading  can  be  related  to  the  average  of  DGD  values. 

In  long  fibre  spans,  DGD  is  random  in  both  time  and  wavelength  since  it  depends  on  the 
details  of  the  birefringence  along  the  entire  fibre  length.  It  is  also  sensitive  to  time-dependent 
temperature  and  mechanical  perturbations  on  the  fibre.  For  this  reason,  a  useful  way  to 
characterize  PMD  in  long  fibres  is  in  terms  of  the  expected  value,  <At>,  or  the  mean  DGD 
over  wavelength.  In  principle,  the  expected  value  <At>  does  not  undergo  large  changes  for  a 
given  fibre  from  day  to  day  or  from  source  to  source,  unlike  the  parameters  5t  or  Ar.  In 
addition,  <A7>  is  a  useful  predictor  of  lightwave  system  performance. 

The  term  "PMD"  is  used  both  in  the  general  sense  of  two  polarization  modes  having  different 
group  velocities,  and  in  the  specific  sense  of  the  expected  value  <At>.  The  DGD  Aror  pulse 
broadening  (Jr  can  be  averaged  over  wavelength,  yielding  <At>;^,  or  time,  yielding  <AT>t,  or 
temperature,  yielding  <AT>-p.  For  most  purposes,  it  is  not  necessary  to  distinguish  between 
these  various  options  for  obtaining  <At>. 

The  coupling  length  is  the  length  of  fibre  or  cable  at  which  appreciable  coupling  between 
the  two  polarization  states  begins  to  occur.  If  the  fibre  length  L  satisfies  the  condition  L  «  1^, 
mode  coupling  is  negligible  and  <At>  scales  with  fibre  length.  The  corresponding  PMD 
coefficient  is 

"short-length"  PMD  coefficient  =  <At>/L. 

Fibres  in  practical  systems  are  nearly  always  in  the  L  »  /^,  regime  and  mode  coupling  is 
random.  If  mode  coupling  is  also  found  to  be  random,  <At>  scales  with  the  square  root  of 
fibre  length,  and 


"long-length"  PMD  coefficient  =  <At>l 4l 
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Indian  Standard 
OPTICAL  FIBRES 

PART  1  MEASUREMENT  METHODS  AND  TEST  PROCEDURES 
Section  48  Polarization  Mode  Dispersion 

1  Scope 

This  part  of  lEC  60793  applies  to  three  methods  of  measuring  polarization  mode  dispersion 
(PMD),  which  are  described  in  Clause  4.  It  establishes  uniform  requirements  for  measuring 
the  PMD  of  single-mode  optical  fibre,  thereby  assisting  in  the  inspection  of  fibres  and  cables 
for  commercial  purposes. 

2  Normative  references 

The  following  referenced  documents  are  indispensable  for  the  application  of  this  document. 
For  dated  references,  only  the  edition  cited  applies.  For  undated  references,  the  latest  edition 
of  the  referenced  document  (including  any  amendments)  applies. 

lEC  60793-1-1 ,  Optical  fibres  -  Part  1-1:  Measurement  methods  and  test  procedures  - 
General  and  guidance 

lEC  60793-1-44,  Optical  fibres  -  Part  1-44:  Measurement  methods  and  test  procedures  - 
Cut-off  wavelength 

I  EC  60793-2-50,  Optical  fibres  -  Part  2-50:  Product  specifications  -  Sectional  specification  for 
class  B  single-mode  fibres 

I  EC  60794-3,  Optical  fibre  cables  -  Part  3:  Sectional  specification  -  Outdoor  cables 

lEC  61280-4-4,  Fibre  optic  communication  subsystem  test  procedures  -  Part  4-4:  Cable 
plants  and  links  -  Polarization  mode  dispersion  measurement  for  installed  links 

lEC/TR  61282-3,  Fibre  optic  communication  system  design  guides  -  Part  3:  Calculation  of  link 
polarization  mode  dispersion 

IEC/TR  61282-9,  Fibre  optic  communication  system  design  guides  -  Part  9:  Guidance  on 
polarization  mode  dispersion  measurements  and  theory 

lEC  61290-11-1,  Optical  amplifier  test  methods  -  Part  11-1:  Polarization  mode  dispersion  - 
Jones  matrix  eigenanalysis  method  (JME) 

lEC  61290-11-2,  Optical  amplifiers  -  Test  methods  -  Part  11-2:  Polarisation  mode  dispersion 
parameter  -  Poincare  sphere  analysis  method 

IEC/TR  61292-5,  Optical  amplifiers  -  Part  5:  Polarization  mode  dispersion  parameter  - 
General  information 

lEC  61300-3-32,  Fibre  optic  interconnecting  devices  and  passive  components  -  Basic  test 
and  measurement  procedures  -  Part  3-32:  Examinations  and  measurements  -  Polarization 
mode  dispersion  measurement  for  passive  optical  components 

ITU-T  Recommendation  G.650.2,  Definitions  and  test  methods  for  statistical  and  non-linear 
related  attributes  of  single-mode  fibre  and  cable 
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3  Terms  and  definitions 

For  the  purposes  of  this  document,  the  terms  and  definitions  contained  in  ITU-T 
Recommendation  G.650.2  apply. 

NOTE    Further  explanation  of  their  use  in  this  document  is  provided  in  lEC  61282-9. 

4  General 

4.1      Methods  for  measuring  PMD 

Three  methods  are  described  for  measuring  PMD  (see  Annexes  A,  B  and  C  for  more  details). 
The  methods  are  listed  below  in  the  order  of  their  introduction.  For  some  methods,  multiple 
approaches  of  analyzing  the  measured  results  are  also  provided. 

-  Method  A 

•  Fixed  analyser  (FA) 

•  Extrema  counting  (EC) 

•  Fourier  transform  (FT) 

•  Cosine  Fourier  transform  (CFT) 

-  Method  B 

•  Stokes  parameter  evaluation  (SPE) 

•  Jones  matrix  eigenanalysis  (JME) 

•  Poincare  sphere  analysis  (PSA) 

•  State  of  polarization  (SOP) 

-  Method  C 

•  Interferometry  (INTY) 

•  Traditional  analysis  (TINTY) 

•  General  analysis  (GINTY) 

The  PMD  value  is  defined  in  terms  of  the  differential  group  delay  (DGD),  At,  which  usually 
varies  randomly  with  wavelength,  and  is  reported  as  one  or  another  statistical  metric. 
Equation  (1)  is  a  linear  average  value  and  is  used  for  the  specification  of  optical  fibre  cable. 
Equation  (2)  is  the  root  mean  square  value  which  is  reported  by  some  methods.  Equation  (3) 
can  be  used  to  convert  one  value  to  the  other  if  the  DGDs  are  assumed  to  follow  a  Maxwell 
random  distribution. 


NOTE  Equation  (3)  applies  only  when  the  distribution  of  DGDs  is  Maxwellian,  for  instance  when  the  fibre  is 
randomly  mode  coupled.  The  generalized  use  of  Equation  (3)  can  be  verified  by  statistical  analysis.  A  Maxwell 
distribution  may  not  be  the  case  if  there  are  point  sources  of  elevated  birefringence  (relative  to  the  rest  of  the 
fibre),  such  as  a  tight  bend,  or  other  phenomena  that  reduce  the  mode  coupling,  such  as  a  continual  reduced  bend 
radius  with  fibre  in  tension.  In  these  cases,  the  distribution  of  the  DGDs  will  begin  to  resemble  the  square  root  of  a 
non-central  Chi-square  distribution  with  three  degrees  of  freedom.  For  these  cases,  the  PMOr^s  value  will 
generally  be  larger  relative  to  the  PMO^vg  ^^^^  indicated  in  Equation  (3).  Time  domain  methods  such  as  Method 
C  and  Method  A,  cosine  Fourier  transform,  which  are  based  on  PMOr^s.  can  use  Equation  (3)  to  convert  to 
PMDavg-  'f  niode  coupling  is  reduced,  the  resultant  reported  PMD  value  from  these  methods  may  exceed  those 
that  can  be  reported  by  the  frequency  domain  measurements  that  report  PMO^vc  such  as  Method  B. 


(1) 


(2) 


(3) 
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The  PMD  coefficient  is  the  PMD  value  normalized  to  the  fibre  length.  For  normal  transmission 
fibre,  for  which  random  mode  coupling  occurs  and  for  which  the  DGDs  are  distributed  as 
Maxwell  random  variables,  the  PMD  value  is  divided  by  the  square  root  of  the  length  and  the 
PMD  coefficient  is  reported  in  units  of  ps/km''^^  por  some  fibres  with  negligible  mode 
coupling,  such  as  polarization  maintaining  fibre,  the  PMD  value  is  divided  by  the  length  and 
the  PMD  coefficient  is  reported  in  units  of  ps/km. 

All  methods  are  suitable  for  laboratory  measurements  of  factory  lengths  of  optical  fibre  and 
optical  fibre  cable.  For  all  methods,  changes  in  the  deployment  of  the  specimen  can  alter  the 
results.  For  installed  lengths  of  optical  fibre  cable  that  may  be  moving  or  vibrating,  either 
Method  C  or  Method  B  (in  an  implementation  capable  of  millisecond  measurement  time 
scales)  is  appropriate. 

All  methods  require  light  sources  that  are  controlled  at  one  or  more  states  of  polarization 
(SOPs).  All  methods  require  injecting  light  across  a  broad  spectral  region  (i.e.  50  nm  to 
200  nm  wide)  to  obtain  a  PMD  value  that  is  characteristic  of  the  region  (i.e.  1  300  nm  or 
1  550  nm).  The  methods  differ  in: 

a)  the  wavelength  characteristics  of  the  source; 

b)  the  physical  characteristics  that  are  actually  measured; 

c)  the  analysis  methods. 

Method  A  measures  PMD  by  measuring  a  response  to  a  change  of  narrowband  light  across  a 
wavelength  range.  At  the  source,  the  light  is  linearly  polarized  at  one  or  more  SOPs.  For  each 
SOP,  the  change  in  output  power  that  is  filtered  through  a  fixed  polarization  analyser,  relative 
to  the  power  detected  without  the  analyser,  is  measured  as  a  function  of  wavelength.  The 
resulting  measured  function  can  be  analysed  in  one  of  three  ways. 

-  By  counting  the  number  of  peaks  and  valleys  (EC)  of  the  curve  and  application  of  a 
formula  that  has  been  shown  [1]'')  to  agree  with  the  average  of  DGD  values,  when  the 
DGDs  are  distributed  as  Maxwellian.  This  analysis  is  considered  as  a  frequency  domain 
approach. 

-  By  taking  the  FT  of  the  measured  function.  This  FT  is  equivalent  to  the  pulse  spreading 
obtained  by  the  broadband  transmission  of  Method  0.  Appropriate  characterisation  of  the 
width  of  the  FT  function  agrees  with  the  average  of  DGD  values,  when  the  DGDs  are 
distributed  as  Maxwellian. 

-  By  taking  the  cosine  Fourier  transform  of  the  difference  of  the  normalized  spectra  from  two 
orthogonal  analyzer  settings  and  calculating  the  RMS  of  the  squared  envelope.  The 
PMD^^^  value  is  reported.  This  is  equivalent  to  simulating  the  fringe  pattern  of  the  cross- 
correlation  function  that  would  result  from  interferometric  measurements. 

Method  B  measures  PMD  by  measuring  a  response  to  a  change  of  narrowband  light  across  a 
wavelength  range.  At  the  source,  the  light  is  linearly  polarized  at  one  or  more  SOPs.  The 
Stokes  vector  of  the  output  light  is  measured  for  each  wavelength.  The  change  of  these 
Stokes  vectors  with  angular  optical  frequency,  ^  and  with  the  (optional)  change  in  input  SOP 
yields  the  DGD  as  a  function  of  wavelength  through  relationships  that  are  based  on  the 
following  definitions: 


6s{cd) 
dco 


Q{co)xs{co) 


(4) 


(5) 


where 


s    is  the  normalized  output  Stokes  vector; 


)  Figures  in  square  brackets  refer  to  the  Bibliography. 
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Q   is  the  polarization  dispersion  vector  (PDV)  in  the  direction  of  the  PSPs; 
^T  is  the  DGD. 

For  both  the  JIVIE  and  PSA  analysis  approaches,  three  linear  SOPs  at  nominally  0°,  45°,  and 
90°  (orthogonal  on  the  Poincare  sphere)  must  be  launched  for  each  wavelength. 

The  JME  approach  is  completed  by  transforming  the  output  Stokes  vectors  to  Jones  matrices 
[2],  appropriate  combination  of  the  matrices  at  adjacent  wavelengths,  and  a  calculation  using 
the  eigenvalues  of  the  result  to  obtain  the  DGD,  by  application  of  an  argument  formula,  at  the 
base  frequency. 

The  PSA  approach  is  completed  by  doing  matrix  algebra  on  the  normalized  output  Stokes 
vectors  to  infer  the  rotation  of  the  output  Stokes  vector  on  the  Poincare  sphere  at  two 
adjacent  wavelengths,  using  the  application  of  an  arcsine  formula  to  obtain  the  DGD.  The 
JME  and  PSA  approaches  are  mathematically  equivalent  for  common  assumptions  (see 
lEC  61282-9). 

The  SOP  approach  is  based  on  a  piecewise  evaluation  of  Equation  (4)  using  the  normalized 
measured  Stokes  vectors.  The  SOP  approach  can  yield  good  results  when  the  transit  of  the 
output  Stokes  vector  is  well  behaved  (negligible  mode-coupling)  but  can  produce  incorrect 
results  when  the  output  Stokes  vector  changes  rapidly  and  randomly  (see  lEC  61282-9).  The 
extra  measurement  time  required  for  the  three  input  SOPs  for  JME  and  PSA  result  in  a  more 
robust  measurement. 

Method  0  is  based  on  a  broadband  light  source  that  is  linearly  polarized.  The  cross- 
correlation  of  the  emerging  electromagnetic  field  is  determined  by  the  interference  pattern  of 
output  light,  i.e.  the  interferogram.  The  determination  of  the  PMD  delay  for  the  wavelength 
range  associated  with  the  source  spectrum  is  based  on  the  envelope  of  the  fringe  pattern  of 
the  interferogram.  Two  analyses  are  available  to  obtain  the  PMD  delay  (see  lEC  61282-9), 
both  of  which  measure  the  PMD^y^^  value: 

-  TINTY  uses  a  set  of  specific  operating  conditions  for  its  successful  applications  and  a 
basic  setup; 

-  GINTY  uses  no  limiting  operating  conditions  but  in  addition  to  the  same  basic  set-up  also 
using  a  modified  setup  compared  to  TINTY. 

With  the  exception  of  the  Method  B  SOP  approach,  the  analysis  approaches  represent  an 
evolution  of  the  understanding  of  PMD.  The  GINTY  is,  for  example,  more  complete  than 
TINTY.  The  reproducibility  of  PMD  depends  on  the  PMD  level  and  the  wavelength  range  of 
the  measurement  [3].  Better  relative  reproducibility  is  achieved  for  broader  wavelength  ranges 
and  higher  PMD  values  for  a  given  range.  For  measurements  of  higher  PMD  values,  e.g., 
0,5  ps,  the  differences  in  the  analysis  methods  are  less  important  than  for  the  measurements 
of  low  PMD  values. 

Information  common  to  all  three  methods  is  contained  in  Clauses  4  to  10,  and  requirements 
pertaining  to  each  individual  method  appear  in  Annexes  A,  B,  and  0,  respectively.  lEC  61282-9 
provides  the  mathematical  formulations  for  all  methods. 

4.2  Reference  test  method 

Method  B,  SPE  (only  JME  and  PSA  approaches),  is  the  reference  test  method  (RTM),  which 
shall  be  the  one  used  to  settle  disputes. 

4.3  Applicability 

PMD  in  fibre  is  a  statistical  parameter.  lEC  60794-3  includes  a  statistical  requirement  on 
PMD,  called  PMDq  or  link  design  value,  that  is  based  on  sampled  measurements  of  optical 
fibre  cable  and  calculations  for  concatenated  links.  The  PMD  of  a  cabled  fibre  can  vary  from 
the  PMD  of  the  uncabled  fibre  due  to  effects  of  cable  construction  and  processing.  A  limit  on 
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the  PMDq  of  the  uncabled  fibre  is,  however,  required  to  limit  the  PMDq  on  cabled  fibre. 
Uncabled  fibre  PMDq  less  than  half  the  cabled  fibre  PMDq  limit  is  generally  considered  as  a 
conservative  rule.  Alternative  limits  may  be  determined  for  particular  constructions  and  stable 
cable  processes. 

The  fibre  or  cable  deployment  should  be  selected  so  externally  induced  mode-coupling  is 
minimized.  Sources  of  such  external  mode-coupling  can  be: 

a)  excessive  tension; 

b)  excessive  bending  induced  from 

•  fibre  cross-overs  on  a  shipping  reel; 

•  crimping  of  fibre  within  a  cable  on  a  spool  that  is  too  small; 

•  too  small  a  bend  radius; 

c)  excessive  twist. 

Reproducibility  of  individual  measurements  should  be  evaluated  after  perturbing  the  fibre  to 
allow  sampling  the  full  range  of  mode-coupling  combinations.  This  can  be  done  by,  for 
example,  changing  the  temperature  slightly  or  making  small  adjustments  in  the  deployment. 
Gisin  [3]  reported  a  fundamental  relative  reproducibility  limit  for  measurements  and  showed 
that  the  relative  reproducibility  improves  as  the  PMD  increases  and  as  the  spectral  width  of 
the  source  increases.  When  PMD  measurements  are  combined  to  evaluate  the  statistical 
specification  of  optical  fibre  cable  (see  lEC  60794-3),  this  variability  leads  to  a  possible 
overstatement  of  the  link  design  value. 

Guidelines  for  the  calculation  of  PMD  for  systems  that  include  other  components  such  as 
dispersion  compensators  or  optical  amplifiers  are  given  in  lEC  61282-3.  Test  methods  for 
optical  amplifiers  are  given  in  lEC  61290-11-1  and  lEC  61290-11-2,  and  other  design  guides 
in  lEC  61292-5.  Test  methods  for  testing  links  including  amplified  ones  are  given  in 
lEC  61280-4-4.  Test  methods  for  optical  components  are  given  in  lEC  61300-3-32.  General 
information  about  PMD,  the  mathematical  formulation  related  to  the  application  of  the  present 
methods,  and  some  considerations  related  to  the  sampling  theory  related  to  the  use  of 
different  light  sources  and  detection  systems  are  given  in  lEC  61282-9. 

5  Apparatus 

The  following  apparatus  is  common  to  all  three  measurement  methods.  Annexes  A,  B,  and  C 
include  layout  drawings  and  other  equipment  requirements  for  each  of  the  three  methods, 
respectively. 

5.1  Light  source  and  polarizers 

See  Annexes  A,  B,  and  C  for  detailed  options  of  the  spectral  characteristics  of  the  light 
source.  The  source  shall  produce  sufficient  radiation  at  the  intended  wavelength(s)  and  be 
stable  in  intensity  over  a  time  period  sufficient  to  perform  the  measurement.  lEC  61282-9 
provides  additional  guides  concerning  the  source  input  SOP,  degree  of  polarization  (OOP), 
use  of  polarizers  and  polarization  controllers. 

5.2  Input  optics 

An  optical  lens  system  or  fibre  pigtail  may  be  employed  to  excite  the  specimen.  It  is  recom- 
mended that  the  power  coupled  into  the  specimen  be  relatively  insensitive  to  the  position  of 
its  input  end  face.  This  can  be  accomplished  by  using  a  launch  beam  that  spatially  and 
angularly  overfills  the  input  end  face. 

If  using  a  butt  splice,  employ  index-matching  material  between  the  fibre  pigtail  and  the 
specimen  to  avoid  interference  effects.  The  coupling  shall  be  stable  for  the  duration  of  the 
measurement. 
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5.3  Input  positioner 

Provide  means  of  positioning  the  input  end  of  the  specimen  to  the  light  source.  Examples 
include  the  use  of  x-y-z  micropositioner  stages,  or  mechanical  coupling  devices  such  as 
connectors,  vacuum  splices,  three-rod  splices,  etc.  The  position  of  the  fibre  shall  remain 
stable  over  the  duration  of  the  measurement. 

5.4  Cladding  mode  stripper 

Use  a  device  that  extracts  cladding  modes.  Under  some  circumstances  the  fibre  coating  will 
perform  this  function. 

5.5  High-order  mode  filter 

Use  a  means  to  remove  high-order  propagating  modes  in  the  desired  wavelength  range  that  is 
greater  than  or  equal  to  the  cut-off  wavelength  (see  lEC  60793-1-44)  of  the  specimen.  For 
example,  a  one-turn  bend  of  radius  =  30  mm  on  the  fibre  is  generally  sufficient. 

5.6  Output  positioner 

Provide  a  suitable  means  for  aligning  the  fibre  output  end  face  to  the  output  optics.  Such 
coupling  may  include  the  use  of  lenses,  or  may  be  a  mechanical  connector  to  a  detector 
pigtail. 

Provide  means  such  as  a  side-viewing  microscope  or  camera  with  a  crosshair  to  locate  the 
fibre  at  a  fixed  distance  from  the  output  optics.  It  may  be  sufficient  to  provide  only  longitudinal 
adjustment  if  the  fibre  is  constrained  in  the  lateral  plane  by  a  device  such  as  a  vacuum  chuck. 

5.7  Output  optics 

See  Annex  A,  B,  or  C,  as  appropriate. 

5.8  Detector 

For  signal  detection,  an  optical  detector  is  used  which  is  linear  and  stable  over  the  range  of 
intensities  and  measurement  times  that  are  encountered  in  performing  the  measurement. 
A  typical  system  might  include  synchronous  detection  by  a  chopper/lock-in  amplifier,  an 
optical  power  meter,  optical  spectrum  analyser,  or  a  polarimeter.  To  use  the  entire  spectral 
range  of  the  source,  the  detection  system  must  have  a  wavelength  range  which  includes  the 
wavelengths  produced  by  the  light  source.  See  Annex  A,  B,  or  C,  as  appropriate,  for 
additional  details. 

5.9  Computer 

Use  a  computer  to  perform  operations  such  as  controlling  the  apparatus,  taking  intensity 
measurements,  and  processing  the  data  to  obtain  the  final  results. 

6    Sampling  and  specimens 
6.1  General 

A  specimen  is  a  known  length  of  single-mode  optical  fibre  (lEC  60793-2-50)  which  may  or 
may  not  be  cabled.  The  sample  and  pigtails  must  be  fixed  in  position  at  a  nominally  constant 
temperature  throughout  the  measurement.  Standard  ambient  conditions  shall  be  employed 
unless  otherwise  specified.  In  the  case  of  installed  fibres  and  cables,  prevailing  deployment 
conditions  may  be  used. 

Mechanical  and  temperature  stability  of  the  test  device  may  be  observed  by  the  following 
procedures.  For  Method  A,  the  output  power  from  the  fibre  at  a  fixed  wavelength  is  measured 
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with  the  output  analyser  in  place.  In  a  time  period  corresponding  to  a  typical  complete 
measurement,  the  output  power  change  should  be  small  relative  to  the  changes  produced  by 
a  wavelength  increment.  For  Method  B,  the  output  SOP  of  the  test  fibre  on  a  Poincare  sphere 
display  is  viewed.  In  a  time  period  corresponding  to  an  adjacent  pair  of  Jones  matrix 
measurements,  the  output  SOP  change  should  be  small  relative  to  the  change  produced  by  a 
wavelength  increment.  Method  0  is  normally  robust  with  regard  to  slight  temperature  change 
or  fibre  movements. 

End  faces  for  the  input  and  output  ends  of  the  test  sample  must  be  prepared  as  appropriate 
for  the  requirements  of  the  apparatus  and  procedure.  Precautions  shall  be  taken  to  avoid  any 
reflections. 

6.2  Specimen  length 

The  specimen  length  is  dictated  by  three  factors: 

a)  minimum  desired  PMD  coefficient; 

b)  mode-coupling  regime; 

c)  signal  to  noise  ratio. 

Each  test  method  and  implementation  is  limited  to  a  minimum  PMD  value  (ps)  that  can  be 
measured.  In  many  cases,  this  minimum  can  be  determined  on  the  basis  of  theory.  It  can  also 
be  determined  experimentally  by  examining  the  measured  distribution.  For  fibres  in  the 
random  mode-coupling  regime,  the  minimum  PMD  coefficient  is  determined  by  dividing  the 
PMD  value  by  the  square  root  of  the  fibre  length  (km).  For  the  negligible  mode-coupling  case, 
the  division  is  by  the  length.  The  length  that  is  measured  and  the  minimum  measurable  PMD 
value  will  therefore  determine  the  minimum  measurable  PMD  coefficient.  Fibres  or  cables  with 
lengths  sufficient  to  achieve  this  minimum  can  be  selected  for  measurement.  Alternatively, 
specimens  can  be  cut  to  a  length  that  is  satisfactory.  The  minimum  measurable  PMD  value 
shall  be  documented.  The  length  of  the  individual  specimens  shall  be  recorded. 

NOTE    The  length  may  also  be  limited  by  the  deployment  method  (see  6.3)  and  instrument  dynamic  range. 

The  values  specified  in  lEC  60794-3  and  lEC  60793-2-50  express  the  PMD  coefficient  in 
terms  of  ps/Vkm  -  in  effect,  these  documents  assume  that  the  length  measured  is  sufficient  to 
induce  the  randomly  mode-coupled  regime.  For  a  given  fibre  type  or  cable  construction,  this 
can  be  confirmed  by  doing  a  cut-back  experiment  in  which  the  PMD  value  is  measured  on  a 
specimen  at  each  of  several  lengths  -  achieved  by  cutting  the  specimen  back  between 
measurements.  Lengths  above  which  there  is  a  square  root  dependence  of  the  PMD  value  on 
length  may  be  considered  as  randomly  mode-coupled. 

The  dynamic  range  is  limited  by  the  method,  the  source  power,  and  the  overall  loss  of  the 
specimen,  which  is  affected  by  length.  This  limit  must  generally  be  determined  on  the  basis  of 
specific  implementations  by  experimental  means.  This  limit  shall  be  documented. 

6.3  Deployment 

The  deployment  of  the  fibre  or  cable  can  influence  the  result.  For  normal  measurements  to  be 
used  in  specification  conformance  evaluation,  the  following  requirements  apply. 

6.3.1      Uncabled  fibre 

It  is  important  to  minimize  deployment  induced  mode  coupling  when  measuring  uncabled 
fibres,  which  is  done  in  order  to  support  the  primary  requirements  of  cabled  fibre  PMDq.  In 
this  case,  the  fibre  shall  be  supported  in  some  manner  (usually  on  a  reel  having  a  minimum 
wind  radius  of  150  mm),  with  essentially  zero  fibre  tension  (typically  less  than  5  g),  and  no 
tensioned  crossovers.  These  deployment  requirements  can  limit  the  length  that  can  be 
measured,  depending  on  the  spool  diameter,  and  can  make  the  measurement  a  destructive 
one.  Multi-layer  windings  are  possible,  but  should  be  qualified  by  comparison  with  single-layer 
results  on  shorter  lengths. 
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The  measurement  of  uncabled  fibre  deployed  on  shipping  spools  is  not  recommended.  PMD 
results  with  this  deployment  have  been  shown  to  be  substantially  less  than  what  would  be 
obtained  in  cable  form  for  high  PMD  fibre  and  substantially  greater  than  what  would  be 
obtained  in  cable  form  for  low  PMD  fibre. 

6.3.2      Optical  fibre  cable 

PMD  measurements  on  fibres  in  cables  wound  on  shipping  drums  may  not  always  reflect  the 
functionally  relevant  PMD  values  for  fibres  in  the  installed  cable  deployment  configuration. 
Consequently,  to  demonstrate  compliance  with  the  cabled-fibre  PMD  specification,  alternative 
deployment  configurations  or  mapping  functions  relating  on-drum  PMD  value  to  off-drum  PMD 
value  may  be  used  for  factory  measurements.  The  exact  deployment  configuration  shall  be 
agreed  upon  between  the  supplier  and  the  customer. 

7  Procedure 

7.1  Deploy  the  fibre  or  cable  and  prepare  the  ends. 

7.2  Attach  the  ends  to  the  input  and  output  optics. 

7.3  Engage  the  computer  to  complete  the  scans  and  measurements  found  in  Annexes  A,  B, 
and  C  for  the  three  measurement  methods. 

7.4  Complete  documentation. 

8  Calculation  or  interpretation  of  results 

Annexes  A,  B,  and  C  provide  calculations  to  convert  the  measured  data  into  PMD  values.  The 
calculation  of  the  PMD  coefficient  is  carried  out  according  to  whether  random  mode  coupling 
or  negligible  mode  coupling  is  present.  For  the  fibres  specified  in  lEC  60793-2-50,  the  PMD 
value  is  normalized  by  the  square  root  of  the  fibre  length  in  units  of  ps/km''^^ 

9  Documentation 

9.1  Information  required  for  each  measurement 

a)  Specimen  identification 

b)  Testing  date 

c)  Specimen  length 

d)  Wavelength  region  (for  example,  1  550  nm) 

e)  PMD  in  units  of  ps,  and  whether  PMD^^g     ^^^rms  reported 

f)  PMD  coefficient  and  its  units  (ps/Vkm  or  ps/km) 

9.2  Information  to  be  available 

a)  Measurement  method  used 

b)  Calculation  approach  used 

c)  Description  of  the  deployment  method  (including  any  fibre  support  mechanism) 

d)  Wavelength  range  used 

e)  For  Methods  A  and  B  with  a  narrowband  source  and  a  step  mode,  the  number  of 
wavelengths  sampled 

f)  For  Method  C,  the  type  of  fringe-detection  technique 
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g)  Description  of  the  equipment 

h)  Date  of  latest  calibration 

i)  Evidence  supporting  the  mode-coupling  regime  (indicated  by  units  of  the  PMD  coefficient) 
j)  For  Method  B  with  narrowband  source  and  a  step  mode,  the  wavelength  range  resolution 
k)  For  Method  B  with  broadband  source  (BBS),  the  centre  wavelength  and  -3  dB  linewidth 

10  Specification  information 

a)  Type  of  fibre  or  cable 

b)  Failure  or  acceptance  criteria 

c)  Wavelength  region 

d)  Any  deviations  from  this  procedure 
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Annex  A 

(normative) 

Fixed  analyser  measurement  method 


This  annex  contains  requirements  specific  to  Method  A  (FA). 
A.1  Apparatus 

Figure  A.1  shows  possible  block  diagrams. 
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Figure  A.1a  -  Narrowband  source 
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Figure  A.1b  -  Broadband  source 
Figure  A.1  -  Block  diagrams  for  IVIethod  A 


A.1.1       Light  source 

In  all  cases,  two  kinds  of  light  sources  may  be  used,  depending  on  the  type  of  analyser. 
A  narrowband  source  such  as  the  broadband  lamp  and  monochromator  combination  shown  in 
Figure  A. 1a  can  be  used  with  a  polarization  analyser.  A  BBS,  shown  in  Figure  A. lb,  can  be 
used  with  a  narrow  bandpass  filtering  analyser,  such  as  an  optical  spectrum  analyser  or  an 
interferometer  used  as  a  FT  spectrum  analyser  placed  before  the  analyser.  In  the  case  of 
BBS,  the  width  of  the  filter  is  taken  as  the  spectral  width  for  the  purpose  of  calculations. 

In  both  cases,  the  spectral  width  shall  be  sufficiently  small  to  maintain  the  desired  degree  of 
polarization  (see  5.1).  In  both  cases,  the  range  of  wavelengths  shall  be  sufficient  to  provide  a 
PMD  measurement  of  sufficient  precision  at  the  specified  wavelength  region  (see  Clause  A. 3). 
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To  insure  that  all  features  in  the  optical  spectrum  are  adequately  resolved,  the  spectral  width 
should  satisfy 

A  Q 

<{8vAt^J-'  (A.1) 

A) 

where 

V  =  c/X    is  the  optical  frequency; 

is  the  spectral  width; 
Ar^gx     is  the  maximum  anticipated  DGD. 

For  A  in  the  vicinity  of  1  550  nm,  Equation  (A.1)  reduces  to  the  condition  that  AA(nm)  should 
be  less  than  the  reciprocal  of  Ar(ps). 

A.1. 2  Analyser 

The  angular  orientation  of  the  analyser  is  not  critical  but  should  remain  fixed  throughout  the 
measurement.  With  negligible  mode-coupling  or  low  PMD  values,  some  adjustment  of  the 
analyser  may  be  helpful  in  maximising  the  amplitude  of  the  oscillations  in  Figure  A. 2  -  which 
can  also  be  achieved  by  rotating  the  fibre  at  splices  or  connectors.  For  the  CFT  approach,  the 
analyzer  must  be  capable  of  being  rotated  to  a  setting  that  is  orthogonal  to  the  initial  setting. 

NOTE    The  analyser  can  be  replaced  by  a  polarimeter. 

A.2  Procedure 


A.2.1      Wavelength  range  and  increment 

The  procedure  requires  measuring  the  power  as  a  function  of  wavelengths  (or  optical 
frequencies)  over  a  range  at  a  defined  wavelength  or  optical  frequency  increment  once  with 
the  analyser  in  the  optical  path  and  once  without  -  or  once  with  the  analyzer  in  the  optical 
path  and  once  with  the  analyzer  set  to  a  position  orthogonal  to  the  initial  setting.  The 
wavelength  range  can  influence  the  precision  of  the  result  (see  Clause  A. 3).  The  wavelength 
increment  should  be  selected  to  satisfy  Equation  (A.1),  with  the  wavelength  increment 
replacing  A/l. 

If  the  FT  or  CFT  approaches  are  used,  the  step  size  should  ideally  be  uniform  in  optical 
frequency  and  the  number  of  steps  should  be  a  power  of  2.  The  monochromator  step-size, 
expressed  in  optical  frequency,  5v,  must  be  a  factor  of  two  smaller  than  the  "oscillation 
frequency"  corresponding  to  the  maximum  DGD  measured.  Because  of  the  large  amount  of 
power  outside  the  second  moment  for  randomly  mode-coupled  fibres,  the  Nyquist  condition 
must  be  at  least  three  times  the  frequency  of  the  second  moment  for  the  maximum  anticipated 
DGD.  That  is: 

7:^>AW  (A.2) 


NOTE  1  If,  from  the  FT,  it  is  evident  that  there  is  significant  energy  near  Armax'  the  measurement  should  be 
repeated  with  a  reduced  increment. 

NOTE  2  The  source  spectral  width  is  generally  equal  to,  or  less  than  the  smallest  wavelength  increment.  For 
example,  for  Ar^ax  -  0.67  ps,  a  monochromator  spectral  width  of  2  nm  at  1  550  nm  {5v  -  249  GHz)  is  typical. 


11 


IS/IEC  60793-1-48  :  2007 


A.2.2      Complete  the  scans 

Complete  the  scan  with  the  analyser  in  the  light  path.  Record  the  received  power  as  Pp^(A). 

Remove  the  analyser  from  the  light  path  and  repeat  the  scan.  Record  the  received  power  as 
PTot(^)- 

Calculate  the  power  ratio,  R(/l)  as  follows.  Figure  A. 2  shows  an  example  of  both  negligibly 
and  randomly  mode-coupled  results. 

RW  =  ^^  (A.3) 
^Tot  j 

An  alternative  procedure  is  to  leave  the  analyser  in  place  on  the  second  scan,  but  rotate  it 
90°.  Record  the  power  as  P^(A).  The  formula  for  the  power  ratio  is  then: 


NOTE  1    The  ratio,  Pa^Pq,  could  also  be  used  when  extrema  counting  is  used. 

NOTE  2  If  a  polarimeter  is  used  as  the  detection  element,  the  normalised  Stokes  parameters  are  measured 
versus  wavelength.  The  three  spectral  functions  (one  per  vector  element)  are  independent  of  received  power  and 
correspond  to  three  independent  power  ratio  functions  that  can  be  analysed  in  the  same  way. 
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Figure  A.2a  -Negligible  mode-coupling 
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Figure  A.2b  -  Random  mode-coupling 
Figure  A.2  -  Typical  results  from  Method  A 
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A.3  Calculations 

There  are  three  approaches  of  calculating  PMD  from  the  R  function  that  is  measured: 

-  extrema  counting; 

-  Fourier  transform; 

-  cosine  Fourier  transform. 

A.3.1       Extrema  counting 

The  function,  R(A),  should  be  obtained  at  equally  spaced  wavelength  intervals  from  a 
minimum  wavelength  of  to  a  maximum  wavelength  of  E  is  the  number  of  extrema  (both 
maxima  and  minima)  within  the  window.  Alternatively,  the  wavelength  range  can  be  redefined 
so  that  and  X2  coincide  with  extrema,  in  which  case  E  is  the  number  of  extrema  (including 
and  X2)  minus  one.  The  formula  for  the  PMD  value,  <At>,  is: 


where  c  is  the  speed  of  light  in  vacuum  and  /c  is  a  mode-coupling  factor  which  equals  1,0  in 
the  absence  of  random  mode-coupling  and  0,82  in  the  limit  of  random  mode-coupling. 

If  a  polarimeter  is  used  as  the  detection  element,  the  average  of  the  values  derived  from  the 
three  normalised  Stokes  parameter  responses  is  taken  as  the  final  value  of  PMD. 

In  the  presence  of  noise,  the  extrema  can  be  difficult  to  determine.  One  solution  is  to  fit  the 
data  to  a  running  polynomial  which  can  be  evaluated  for  extrema  at  every  point.  A  cubic 
polynomial  that  covers  eight  wavelengths  has  been  used  successfully. 

A.3. 2      Fourier  transform 

In  this  method  a  Fourier  analysis  of  R(k),  usually  expressed  in  the  optical  frequency  domain  v, 
is  used  to  derive  PMD.  The  Fourier  transform  transforms  this  optical  frequency  domain  data  to 
the  time  domain.  The  Fourier  transform  yields  direct  information  on  the  distribution  of  light 
arrival  times  St.  This  data  is  post-processed  as  described  below  to  derive  the  expected  PMD, 
<At>,  for  the  fibre  under  test.  This  method  is  applicable  to  fibres  with  negligible  or  random 
mode  coupling. 

A.3. 2.1       Data  pre-processing  et  Fourier  transformation 

To  use  this  method,  the  Fourier  transform  normally  requires  equal  intervals  in  optical 
frequency  so  that  R{X)  data  are  collected  at  X  values  such  that  they  form  equal  intervals  in  the 
optical  frequency  domain.  Alternatively,  data  taken  at  equal  X  intervals  may  be  fitted  (for 
example,  by  using  a  cubic  spline  fit)  and  interpolation  used  to  generate  these  points,  or  more 
advanced  spectral  estimation  techniques  can  be  used.  In  each  instance,  the  ratio  R{X)  at  each 
/lvalue  used  is  calculated  using  Equation  (A.3)  or  Equation  (A. 4)  as  appropriate. 

Zero-padding  or  data  interpolation  and  DC  level  removal  may  be  performed  on  the  ratio  data, 
R{X).  Windowing  the  data  may  also  be  used  as  a  pre-conditioning  step  before  the  Fourier 
transform.  The  Fourier  transform  is  now  carried  out,  to  yield  the  amplitude  data  distribution 
P{dT)  for  each  value  of  St. 


kEX  1 X  2 


(A5) 
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A. 3. 2. 2       Transform  data  fitting 

Fourier  transform  data  at  zero  dr  has  little  meaning  since,  unless  carefully  removed,  DC 
components  in  R{X)  may  be  partially  due  to  insertion  loss  of  the  analyser  for  example.  When 
the  DC  level  is  not  removed,  up  to  two  data  points  are  generally  bypassed  (not  used)  in  any 
further  calculations.  A  variable,  y,  can  be  defined  so  that  the  'first  valid  bin'  above  zero  (Jrthat 
is  included  in  calculations  corresponds  to  j  =  0. 

In  order  to  remove  measurement  noise  from  subsequent  calculations,  P{St)  is  compared  to  a 
threshold  level  typically  set  to  200%  of  the  RMS  noise  level  of  the  detection  system.  It  is 
now  necessary  to  determine  whether  the  fibre  is  negligibly  or  randomly  mode  coupled. 

If  it  is  found  that  the  first  X  valid  points  of  P(St)  are  all  below  7^,  this  indicates  that  P(St)  must 
have  discrete  spike  features  characteristic  of  negligibly  coupled  fibres.  The  value  of  X  is  equal 
to  three,  unless  zero-padding  is  used  in  the  Fourier  analysis.  In  that  case,  the  value  of  X  can 
be  determined  from 

3  X  (the  number  of  original  data  points) 
(total  length  of  array  after  zero-padding) 

PMD  is  calculated  using  Equation  (A. 6)  for  a  negligible  mode  coupling  fibre,  or  PMD  is 
calculated  using  Equation  (A. 7)  for  a  random  mode  coupling  fibre. 

A. 3. 2. 2.1        PIVID  calculation  for  fibres  with  negligible  mode  coupling 

For  a  negligibly  coupled  fibre  (e.g.,  a  high  birefringence  fibre)  or  for  a  birefringent  component, 
R(A)  resembles  a  chirped  sine  wave  (Figure  A. 2a).  Fourier  transform  will  give  a  P(St)  output 
containing  a  discrete  spike  at  a  position  corresponding  to  the  relative  pulse  arrival  time,  St, 
the  centroid  of  which  is  the  PMD  value  <At>. 

To  define  the  spike  centroid  <Af>,  those  points  where  P{St)  exceeds  a  second  pre-determined 
threshold  level  T2,  typically  set  to  200  %  of  the  RMS  noise  level  of  the  detection  system,  are 
used  in  the  equation: 

<Ar>=X^^o[^-(^^)^^-l^Sr^o[^-(^^)  (A.6) 

where  M+1  is  the  number  of  data  points  of  P  within  the  spike  which  exceed  T2.  <At>  in 
Equation  (A.6)  is  typically  quoted  in  picoseconds.  If  no  spike  is  detected  (i.e.,  M'  =  0),  then 
PMD  is  zero.  Other  parameters  such  as  the  RMS  spike  width  and/or  spike  peak  value  may  be 
reported. 

If  the  device  under  test  contains  one  or  more  birefringent  elements,  more  than  one  spike  will 
be  generated.  For  a  number  n  concatenated  fibres/devices,  up  to  2^'^"'')  spikes  will  be 
obtained. 

A.3.2.2.2       PMD  calculation  for  fibres  with  random  mode  coupling 

In  instances  of  random  mode  coupling,  R(A)  becomes  a  complex  waveform  similar  to 
Figure  A. 2b,  the  exact  characteristics  being  based  on  the  actual  statistics  of  the  coupling 
process  within  the  fibre/cable.  The  Fourier  transformed  data  now  becomes  a  distribution  P(St) 
representing  the  combination  of  autocorrelation  and  cross-correlation  functions  of  light  pulse 
arrival  times,  St,  in  the  fibre  (see  Figure  A. 3). 
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Counting  up  from  j  =  0,  the  first  point  of  P  is  determined  which  exceeds  7^,  and  which  is 
followed  by  at  least  X data  points  which  fall  below  T^.  This  point  represents  the  last  significant 
point  in  (i.e.,  the  'end'  of)  the  distribution  P{St),  for  a  randomly  mode-coupled  fibre,  that  is  not 
substantially  affected  by  measurement  noise.  The  (Jr  value  for  this  point  is  denoted  Sr^^^^^,  and 
the  value  of  j  at  Sr^^^^^  is  denoted  "M". 

The  square  root  of  the  second  moment,  Op,  of  this  distribution  defines  the  fibre  PMD  <At>, 
and  is  given  by: 

1 

<  At  >=  ctr  =  {Z'^"qIPj{St)St^]I^'^"^[Pj{St)]}'^  (A.7) 
A.3.2.2.3       PMD  calculation  for  mixed  coupling  fibre  systems 

There  may  be  instances  where  both  negligibly  coupled  fibre/components  and  randomly 
coupled  fibre(s)  are  concatenated  to  form  the  system  under  test.  In  this  case,  both  centroid 
determination  (Equation  (A.6))  and  the  second  moment  derivation  (Equation  (A.7))  may  be 
required.  Note  that  spikes  in  P(St)  may  only  be  determined  beyond  the  ^^Tjgg.  computed. 
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Figure  A.3  -  PIVID  by  Fourier  analysis 


A.3.3      Cosine  Fourier  analysis 

This  analysis  is  based  on  the  observation  that  the  cosine  Fourier  transform  of  the  spectrum 
emitted  from  the  analyzer  is  the  fringe  pattern  of  the  interferogram  that  would  be  obtained 
from  Method  C.  The  difference  between  fringe  patterns  generated  by  the  analyzer  being  set  at 
two  orthogonal  settings  is  the  cross-correlation  function.  For  an  infinite  spectrum  into  the 
analyzer,  the  autocorrelation  function  would  have  zero  width.  In  practice,  the  finite  source 
spectrum  in  the  optical  frequency  domain  acts  as  a  windowing  function  which  produces  a  non- 
zero autocorrelation  function  width  in  the  time  domain. 


16 


IS/IEC  60793-1-48:2007 


The  analysis  of  the  squared  cross-correlation  and  autocorrelation  functions  found  in  the 
Method  C,  GINTY  analysis  [4]  shows  that  the  difference  in  squared  RMS  widths  of  these 
functions  is  proportional  to  the  square  of  the  spectrally  weighted  RMS  (by  squared  power)  of 
the  DGD  values.  (See  Equation  (C.9).) 

The  result  is  independent  of  the  spectral  shape  which  means  that  the  details  of  the  windowing 
function  are  fully  taken  into  account.  It  is  also  independent  of  the  degree  of  nnode  coupling, 
which  means  that  no  changes  in  algorithm  are  needed  to  treat  the  different  mode  coupling 
regimes. 

The  result  is  limited  by  the  spectral  width  and  optical  frequency  increment  that  is  measured. 
As  the  PMD  increases,  the  frequency  increment  must  be  decreased.  At  some  limit  it  would  be 
more  practical  to  use  the  Method  C  (GINTY). 

The  analysis  reports  the  PMD^^^  metric.  If  random  mode  coupling  is  present,  the  result  can 
be  converted  to  PMDj^^q  using  Equation  (3). 

A.3.3.1  Overview 

The  measurement  of  the  powers  emitted  from  the  analyzer  set  at  two  orthogonal  settings  is 
required.  The  ratio,  R,  associated  with  Equation  (A.4)  is  modified  to: 

where  v  =  c/ A  is  the  optical  frequency  (THz). 

If  a  polarimeter  is  used,  the  three  normalized  output  Stokes  vector  elements  are  equivalent  to 
three  independent  normalized  ratios  equivalent  to  that  represented  by  Equation  (A. 8).  Each 
Stokes  vector  element  is  the  difference  in  powers  between  orthogonal  analyzer  settings.  The 
three  elements  are  different  in  that  the  base  settings  are  also  orthogonal. 

The  data  is  multiplied  by  a  windowing  function,  W{v),  that  goes  to  zero  smoothly  at  the  edges. 
Both  R(v)W(v)  and  W(v)  are  put  into  arrays  with  zero  padding  at  lower,  unmeasured 
frequencies.  Fast  cosine  Fourier  transforms  (FCFT)  are  applied  to  each  array  to  obtain  the 
time  domain  fringe  envelopes,  r{t)w{t)  and  w{t).  These  are  squared  to  obtain  the  squared 

cross-correlation  and  autocorrelation  envelopes,  ExandEg  ,  respectively.  When  multiple  ratio 

functions  (A/)  are  available  from  different  combinations  of  input  polarizer  setting  and  base 
analyzer  settings  (or  different  Stokes  output  vector  elements),  using  for  instance  input/output 
SOP  scrambling,  form  the  mean  square  envelopes  as: 


^x=^E^x/  (A.9) 

/ 


Using  the  RMS  calculation  of  Clause  D.2,  calculate  the  RMS  widths,  and  Oq  of  these  two 
functions.  The  PMD^^^  value  is  calculated  as: 


PMDi 


RMS 


1/2 

(A.11) 
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It  is  related  to  the  spectrally  weighted  (by  squared  window  value)  RMS  of  the  DGDs  as: 

The  expected  value  operator  is  with  respect  to  random  input/output  SOPs. 
A.3.3.2  Details 

This  subclause  explains  some  of  the  details  with  respect  to  the  measured  frequency  window, 
the  frequency  increment,  Av,  frequency  shifting,  and  the  result  of  the  FCFT.  An  example  of  a 
FCFT  algorithm  may  be  found  in  [5]. 

The  data  must  be  available  in  uniform  frequency  increments.  The  number  of  data  points, 
including  zero  pad  values,  must  be  1+2^,  with  k  an  integer. 

If  the  nm  measured  data  points  are  not  taken  in  uniform  frequency  increments,  they  may  be 
fitted  to  a  polynomial  such  as  a  spline  for  interpolation.  A  cubic  spline  [6]  with  nm  -  3  uniform 
segments  will  fit  all  the  data  perfectly  and  allow  interpolation. 

Given  that  the  measured  data  are  bounded  by  v^jp/u  and  v^gxM  ^^^^  ^^^^  minimum 

optical  frequency  is  well  above  zero,  the  application  of  frequency  shifting  can  be  used  to 
reduce  the  size  of  the  arrays  that  are  processed.  The  boundaries  of  the  frequencies  used  in 
the  calculation  array  can  be  selected  by  any  choice  of  n  such  that: 

A?  -1 

^max— ^  =  ^min  ^^minM-  ^max  =^maxM>  and  A?  is  a  positive  integer.  (A.  13) 

The  frequency  values  less  than  the  measured  frequency  are  filled  with  zeros. 

Following  the  FCFT,  the  array  will  contain  the  time  domain  fringe  pattern  from  times  of  0  to 
fmax  =  ^^2^ ,  where  the  time  increment.  At,  is  given  as: 


A^  =  -  =  -j  Y  (A.14) 

The  fringe  pattern  that  would  be  obtained  from  interferometry  extends  to  negative  time  values 
as  well  as  positive  time  values.  The  value  at  a  given  negative  time  is  equal  to  the  value  at  the 
positive  time.  The  function  is  even  and  symmetric  about  zero. 

The  selection  of  the  frequency  shift  should  be  done  keeping  in  mind  that  the  RMS  width 
calculation  needs  some  time  domain  values  that  are  less  than  the  minimum  PMD^^^^  that  is 
measurable. 

The  frequency  increment,  Av,  is  also  related  to  the  number  of  points  sampled,  the  frequency 
shift,  and  the  maximum  PMD^^^  that  is  to  be  measured.  It  is  given  as  the  following,  along 
with  the  constraint  as: 

Av  =  ^nriax  -^min  ^  1   ^ 

2^  24PMDRMS-max 
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The  spectral  width  of  the  filtered  source  should  be  half  of  this  value.  When  the  actual  scan  is 
done  in  equal  wavelength  increments,  the  wavelength  increment  at  the  lower  end  of  the  range 
should  be  consistent  with  the  constraint  of  Equation  (A. 14). 

The  windowing  function,  W(v),  can  technically  be  any  function,  including  a  square  function. 
The  function  that  is  chosen  should  be  one  that  minimizes  the  value  of  a^.  Functions  that  do 
this  proceed  to  zero  at  the  edges  in  a  continuous  way  and  should  also  have  the  first  derivative 
proceed  to  zero  at  the  edges.  This  will  minimize  the  ringing  that  can  increase  Oq- 

A.3.3.3  Examples 

Table  A.1  shows  a  sample  calculation  spreadsheet.  The  wavelength  extrema  and  k  are 
entered.  For  each  of  several  possible  frequency  shift  values,  n,  the  other  parameters  are 
calculated.  minPMD  is  calculated  as  3Af.  The  increment  in  terms  of  AA  at  the  lower 
wavelength  limit  is  also  presented.  Clearly  there  are  tradeoffs,  depending  on  the  range  of 
PMD^^^  values  that  are  to  be  measured.  In  general,  the  broader  the  wavelength  range  and 
the  smaller  the  frequency  increment,  the  better. 


Table  A.1  -  Cosine  transform  calculations 


c 

299792.5 

nm/ps 

k 

11 

npt 

2048 

wavelengths  (nm) 

frequencies  (THz) 

max 

1700 

230.6096 

min 

1300 

176.3485 

Frequency  shift 

wave 

freqmincalc 

delfreq 

minPMD 

maxP/WD 

increment 

n 

(THz) 

del-t  (ps) 

(THz) 

(ps) 

(ps) 

(nm) 

1 

0 

0.002168 

0.112602 

0.006504 

0.370034 

0.635076 

2 

115.3048077 

0.004336 

0.225205 

0.013009 

0.185017 

1 .270772 

3 

153.7397436 

0.006504 

0.337807 

0.019513 

0.123345 

1 .90709 

4 

172.9572115 

0.008673 

0.450409 

0.026018 

0.092508 

2.544031 

5 

184.4876923 

0.010841 

0.563012 

0.032522 

0.074007 

3.181596 

6 

192.1746795 

0.013009 

0.675614 

0.039027 

0.061672 

3.819785 

7 

197.6653846 

0.015177 

0.788216 

0.045531 

0.052862 

4.458599 

8 

201.7834135 

0.017345 

0.900819 

0.052036 

0.046254 

5.098039 

Figure  A. 4  illustrates  the  results  for  what  might  be  obtained  from  a  fibre  with  PMD^^^  = 
0,2  ps.  The  mean  cross-correlation  and  mean  autocorrelation  envelopes  from  a  single  scan 
using  a  windowing  function  that  is  Gaussian  with  a  standard  deviation  of  23  nm  are  shown. 
The  result  is  from  a  simulation  of  a  fibre  with  ideal  random  mode  coupling.  The  measured 
result  for  this  simulation  was  0,185  ps. 
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Figure  A.4  -  Cross-correlation  and  autocorrelation  functions 
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Annex  B 

(normative) 

Stokes  evaluation  method 


This  annex  contains  requirements  specific  to  IVIethod  B  (SPE). 
B.1  Apparatus 

Figure  B.1  shows  a  possible  block  diagram  for  Method  B. 
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Figure  B.1a  -  Narrowband  source  (JME) 
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Figure  B.1b  -  Broadband  source  (PSA) 
Figure  B.1  -  Block  diagram  for  IVIethod  B 


B.1.1       Light  source 


In  all  cases,  two  kinds  of  light  sources  may  be  used,  depending  on  the  type  of  polarimeter. 
A  narrowband  source  such  a  tuneable  laser  shown  in  Figure  B.1a  can  be  used  with  a 
polarization  analyser.  Alternatively,  a  BBS  shown  in  Figure  B.1b  can  be  used  with  a  narrow 
bandpass  filtering  polarimeter  such  as  an  optical  spectrum  analyser  or  an  interferometer  used 
as  a  FT  spectrum  analyser  placed  before  the  polarimeter.  In  the  case  of  BBS,  the  width  of  the 
filter  is  taken  as  the  spectral  width  for  the  purpose  of  calculations. 
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In  both  cases,  the  spectral  width  shall  be  sufficiently  small  to  maintain  the  desired  DOP  (see 
5.1).  In  both  cases,  the  range  of  wavelengths  shall  be  sufficient  to  provide  a  PMD 
measurement  of  sufficient  precision  at  the  specified  wavelength  region  (see  Clause  B.3). 

For  the  JME  and  PSA  approaches,  the  polarizer  must  be  capable  of  switching  between  three 
linear  SOPs  that  are  orthogonal  (nominally  0°,  45°,  and  90°)  for  each  wavelength  measured. 

B.1.2  Polarimeter 

Use  a  polarimeter  to  measure  the  output  Stokes  vectors  for  each  selected  input  SOP  and 
wavelength. 


B.2  Procedure 


The  output  of  the  fibre  is  coupled  to  the  polarimeter.  The  wavelengths  are  scanned  across  a 
range  appropriate  for  the  wavelength  region  and  desired  precision  (see  Clause  B.3)  with  a 
wavelength  increment,  5X.  For  narrowband  sources,  the  wavelength  increment  is  given  in 
terms  of  the  maximum  anticipated  DGD  value,  Ar^g^'  ^^e  wavelength  of  the  region  measured, 
Aq,  and  the  speed  of  light  in  vacuum,  c,  as: 


5X< 


max 


(B.1) 


For  example,  the  product  of  maximum  DGD  and  step  size  shall  remain  less  than  4  ps-nm  at 
1  550  nm  and  less  than  2,8  ps-nm  at  1  300  nm.  This  requirement  ensures  that  from  one  test 
wavelength  to  the  next,  the  output  SOP  rotates  less  than  180°  about  the  PSP  axis  of  the 
Poincare  sphere.  If  a  rough  estimate  of  Ar^^gx  cannot  be  made,  a  series  of  sample 
measurements  is  performed  across  the  wavelength  range,  each  measurement  using  a  closely 
spaced  pair  of  wavelengths  appropriate  to  the  spectral  width  and  minimum  tuning  step  of  the 
optical  source.  The  maximum  DGD  measured  in  this  way  is  multiplied  by  a  safety  factor  of 
three,  this  value  is  substituted  for  Ar^g^  in  the  above  expression  and  the  value  of  5X  to  be 
used  in  the  actual  measurement  is  computed.  If  there  is  concern  that  the  wavelength  interval 
used  for  a  measurement  was  too  large,  the  measurement  may  be  repeated  with  smaller 
wavelength  interval.  If  the  shape  of  the  curve  of  DGD  versus  wavelength  and  the  mean  DGD 
are  essentially  unchanged,  the  original  wavelength  interval  was  satisfactory. 

For  BBSs,  the  resolution  bandwidth  (RBW)  of  the  analyser  must  satisfy  the  following: 


RB\N  <  ^   (B.2) 

The  measurement  data  is  gathered  for  each  wavelength.  For  the  JME  and  PSA  calculation 
approaches,  the  three  input  SOPs  are  cycled  between  nominal  0°,  45°,  and  90°  for  each 
wavelength  are  recorded  in  corresponding  vectors.  For  the  PSA  and  JME  approaches,  the 

output  vectors  are  normalized  to  unit  length  and  recorded  as  H  ,  Q,  and  y  for  the  three  input 
SOPs,  respectively.  For  the  SOP  method,  the  normalized  output  Stokes  vector  for  each 
wavelength  is  recorded  ass. 
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B.3  Calculations 

All  three  calculation  approaches  require  evaluation  of  differences  between  the  SOP  at  one 
angular  optical  frequency,  cdq,  and  the  next  at  cdq  +  Acd.  (Angular  optical  frequency  is  given 
by  CO  =  2nc/ A.)  All  three  calculation  approaches  result  in  a  series  of  DGD  values  versus 

wavelength.  For  negligible  mode-coupling,  the  DGD  values  are  typically  constant  versus 
wavelength.  For  random  mode-coupling,  the  DGD  values  typically  vary  versus  wavelength  as 
shown  in  Figure  B2.  Alternatively,  the  DGD  values  may  be  displayed  as  a  histogram  such  as 
Figure  B.3.  The  average  of  these  DGD  values  is  reported  as  the  PMD  value  that  is  used  in 
conjunction  coefficient. 

The  detailed  mathematical  formulation  and  the  calculations  pertaining  to  Method  B  for  the 
three  approaches,  as  well  as  the  theoretical  linkage  between  the  JME  and  PSA  calculation 
approaches  are  given  in  lEC  61282-9.  The  calculation  for  the  SOP  approach  is  related,  but 
not  as  robust. 

0,2  1  

0,18  


0,16 


0,02 


0  ^  ^  1  ^  \  ^  \  ^  1  ^  1  ^  \  ^  \  ^  \  ^  1 

1  470         1  480         1  490         1  500         1  510         1  520         1  530         1  540         1  550         1  560 
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Figure  B.2  -  Typical  random-mode-coupling  results  from  Method  B 


NOTE    A  Maxwell  curve  is  superimposed  on  the  histogram. 

Figure  B.3  -  Typical  histogram  of  DGD  values 
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The  JME  and  PSA  analysis  approaches  are  mathematically  equivalent  for  first  order 
assumptions  in  the  absence  of  PDL. 

B.3.1      Jones  matrix  eigenanalysis  (JME) 

For  each  frequency  pair,  indicated  by  frequencies,  cdq  and  coq+Aco,  the  three  output  Stokes 
vectors  for  each  frequency  are  converted  to  Jones  vectors  and  a  T  matrix  is  calculated  for 
each  frequency  using  ratios  of  the  elements  of  the  Jones  vectors.  The  following  relationship  is 

used  to  convert  a  normalized  output  Stokes  vector,  noted  generically  as  s,  to  a  Jones  vector, 
noted  generically  as  j  . 


cos  20 
sin2<9cos// 
sin  2(9  sin// 


cos<9exp(-////2) 
sin<9exp(////2) 


(B.3) 


where  0  \s  the  linear  polarization  parameter  and  ju  is  the  circular  parameter,  which  is  also  the 
phase  separation  of  the  x  and  y  element  of  the  Jones  vector.  The  linear  parameter  can  be 
assumed  to  be  within  the  range  0  to  tt  for  this  calculation. 

For  each  frequency,  the  x  and  y  elements  of  the  Jones  vectors  are  designated  as:  h^,  hy,  q^, 
Qy,  v^,  and  Vy.  Using  these,  calculate  the  following  ratios: 


/C-l     =   hx   I  hy  =Vx   I  Vy 


-/C3 


(B.4) 


For  each  frequency,  the  Jones  transfer  matrix,  T,  is  calculated  as: 


7  = 


/C1/C4  /C2 
/C4  1 


(B.5) 


Using  the  T  matrices  for  the  frequency  pair,  the  frequency  transfer  matrix,  J,  is  calculated  as: 

J  =  T{cdq  +Acd)T~\coq)  (B.6) 

The  eigenvalues  of  J  are  calculated  and  designated  as  and  The  DGD,  Ar,  is  calculated 
as: 


Ar  = 


Arg 


Aco 


(B.7) 


where  Arg  denotes  the  argument  function,  that  is,  arg{me^^)=  0 ,  m  >  0  and  \0\  <  k. 


The  JME  mathematical  formulation  and  detailed  calculations  are  given  in  lEC  61282-9. 
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B.3.2      Poincare  sphere  analysis  (PSA) 


For  the  PSA,  matrix  algebra  is  done  on  the  normalized  output  Stokes  vectors  to  deduce  the 
rotation  of  the  output  Stokes  vector  with  frequency. 

For  each  frequency  pair,  indicated  by  frequencies,  co^^  and  co^^+Aco,  the  three  measured  output 
Stokes  vectors  for  each  frequency  are  converted  as  follows: 


h  =  H 


HxQ 


HxQ 


qxV 

v=— — —xq 

qxV 


(B.8) 


From  the  Stokes  vectors  h,  v  and  q  ,  form  the  vector  products  c  =  hxq  and  c'  =  qxv  at  each 
frequency. 

For  each  frequency  pair,  compute  the  finite  differences. 

Ah  =  h(coQ  +  Acd)  -  h(cDQ )       Aq  =  q(coQ  +  Aco)  -  qicoQ )       Av  =  v(coq  +  Aco)  -  v(cdq  ) 

Ac  =  c{cdq  +  Acd)-c{cdq)       Ac'=  c'(^o  +  A^)-c'(^o)  (^-9) 
Find  the  DGD,  Ar ,  for  a  particular  frequency  increment  from  the  following: 


At-- 


1 

Acd 


arcsin 


+  arcsin 


Aq^  +Av^  +Ac'^ 


(B.10) 


where  Ah  =  A/?* A/? 

The  PSA  mathematical  formulation  and  detailed  calculations  are  given  in  lEC  61282-9. 
B.3.3      State  of  polarization  (SOP) 

For  the  SOP  analysis,  the  trace  on  the  Poincare  sphere  describing  the  evolution  of  the  SOP 
with  wavelength  is  reconstructed  from  the  measured  normalized  output  Stokes  vectors.  The 
trace  is  analysed  piecewise,  considering  wavelength  intervals  (which  may  include  more  than 
two  wavelength  steps)  such  that  the  assumptions  ensuring  the  existence  of  well  determined 
PSPs  hold.  The  local  PSP  axis  on  the  Poincare  sphere  and  the  corresponding  rotation  angle 
A<9  caused  by  the  considered  wavelength  variation  SA  are  then  determined  by  means  of  simple 
geometrical  considerations. 

A  possible  procedure  could  be  the  analysis  of  the  trace  on  the  Poincare  sphere  by 
considering  the  measured  points  three  by  three  and  finding  the  point  of  intersection  of  the 
axes  of  the  segments  identified  by  the  two  pairs  of  points.  Starting  from  this  point  it  is 
possible  to  calculate  the  value  of  AO  by  means  of  trigonometric  relationships. 


The  DGD  is  found  from  the  following  expression: 


Af- 


AO 


Aco 


AO 


2nc  ■  SA 


(B.11) 


where  /Ijand      are  respectively  the  initial  and  the  final  wavelength  of  SA. 

The  SOP  mathematical  formulation  and  detailed  calculations  are  given  in  lEC  61282-9. 

NOTE  If  the  output  Stokes  vector  is  aligned  with  the  true  PSP  at  a  given  frequency,  the  calculated  DGD  at  that 
frequency  can  be  substantially  less  than  the  actual  value. 
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Annex  C 

(normative) 

Interferometry  method 


This  annex  contains  detailed  requirements  for  completing  PMD  measurements  using 
Method  C  (INTY). 

C.1  Apparatus 

Figure  C.1  shows  a  block  diagram  for  a  generic  implementation  of  Method  C. 


Optical  frequency  (v)  Delay  (r) 


Figure  C.1  -  Schematic  diagram  for  IVIethod  C  (generic  implementation) 

Parameters  used  in  Figure  C.1  and  later  throughout  the  text: 


V  optical  frequency  (Av  =  c); 

T  difference  of  round-trip  delay  between  the  two  arms  of  the  interferometer; 

Ss(^)  optical  spectrum,  at  FUT  input  =  spectral  density  of  Es(v)the  source  electric  field 
spectrum; 

So(^)  optical  spectrum,  at  FUT  output  (analyser  input); 

S{v)  optical  spectrum,  at  analyser  output  (interferometer  input); 

Sq  input  SOP  (at  FUT  input;  a  unit  Stokes  vector); 

s(v)  output  SOP  (at  FUT  output); 

Sg  analyser  transmission  axis; 
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s(v)*Sa      Stokes  parameter  giving  the  projection  of  s(v)on  the  analyser  transmission  axis. 
It  is  this  parameter  that  contains  the  PMD  information; 

^(^^  optical  power  at  the  interferometer  output,  as  a  function  of  delay  r, 

^(^)  dependent  part  of  P(r)  ("a.c."  part) 

Pq  constant  part  of  P{t)  ("d.c."  parts) 

E(t)  fringe  envelope 

Ex(r)  cross-correlation  envelope 

Eo(x)  autocorrelation  envelope 

The  optical  power  at  the  interferometer  output,  P(t),  is  equal  to  the  sum  of  "a.c."  and  "d.c." 
parts.  Both  parts  are  equal  at  r  =  0  so  the  "a.c."  part  can  be  calculated.  For  an  ideal 
interferometer,  the  "a.c."  part  is  an  even  function,  the  right  half  of  which  is  equal  to  the  cosine 
Fourier  transform  of  the  optical  spectrum,  S(v)^  emitted  from  the  analyser.  For  non-ideal 
interferometers,  some  corrections  may  be  applied,  depending  on  the  details  of  the 
implementation. 


For  TINTY,  the  envelope  of  the  interferogram,  E(t),  is  the  absolute  value  of  the  "a.c."  part.  For 
GINTY,  additional  calculations  to  obtain  the  cross-correlation  and  autocorrelation  envelopes 
are  described  in  C.2.2.2  and  C.3.2.  These  calculations  involve  two  measured  interferograms 
resulting  from  the  analyzer  being  set  at  two  orthogonal  SOPs. 


Figure  C.2  shows  block  diagrams  for  three  specific  implementations. 
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Fringe  envelope 
detection 


Optical 
source 


Fibre 


Polarizer 


Controller 


Connector  or  splice  Mirror  Coupler 


Moving  mirror 


Optical 
source 


AI2 


Polarizer 


Moving  arm 


Controller 


Figure  C.2a  -  Michelson  interferometer 


Beam 
splitter 


Moveable  cube  '^'^I'e 
corner 


Detection 


Connector  or  splice 
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Figure  C.2b  -  IVIach-Zehnder  interferometer 


Polarizer  Scramblers 


DUT 


Mirror 


Beam 
splitter 


Mirror 


Interferometer 


7^ 


Polarization 

diversity 

detector 


Polarization 
beam  splitter 


Fringes:  Px(t)  Pv(t) 


NOTE  The  optional  input/output  (l/0)-SOP  scramblers  are  used  to  improve  the  measurement  uncertainty  (see 
lEC  61282-9). 

Figure  C.2c  -  Setup  with  polarization  scramblers 
Figure  C.2  -  Other  schematic  diagrams  for  IVIethod  C 
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C.1.1       Light  source 

A  BBS  is  used  that  emits  radiation  at  the  intended  nneasurement  wavelengths,  such  as  a  light 
emitting  diode  (LED),  an  amplified  spontaneous  emission  (ASE)  source  or  a  superfluorescent 
source.  The  light  shall  be  polarized  as  shown  in  Figure  C.1.  The  central  wavelength,  Aq,  shall 
be  within  the  1  310  nm  or  1  550  nm  windows  or  any  other  window  of  interest.  In  order  to 
successfully  use  a  TINTY-based  measurement  system,  the  BBS  spectral  shape  shall  be 
approximately  Gaussian,  without  ripples  that  could  influence  the  autocorrelation  function  of 
the  emerging  light.  A  GINTY-based  measurement  system  does  not  require  such  source 
characteristics:  any  shape  can  be  used.  The  spectral  source  line  width  (-3  dB),  AA,  must  be 
known  to  calculate  the  coherence  time,  t^,  which  is  determined  with  the  following: 

AA-  c 


C.1.2      Beam  splitter 

The  beam  splitter  is  used  to  split  the  incident  polarized  light  into  two  components  propagating 
in  the  arms  of  the  interferometer.  The  splitter  can  be  an  optical  fibre  coupler  or  a  cube  beam 
splitter. 

C.1.3  Analyzer 

The  analyzer  function  shown  in  Figure  C.1  may  be  implemented  within  the  interferometer.  For 
the  TINTY  approach,  the  analyzer  must  be  capable  of  being  rotated  to  a  second  setting  that  is 
orthogonal  to  the  initial  setting. 

C.1 .4  Interferometer 

The  interferometer  can  be  an  air  type  or  a  fibre  type.  It  can  be  of  Michelson  or  Mach-Zehnder 
types,  and  it  can  be  located  at  the  source  or  at  the  detector  end  of  the  fibre  under  test.  In  all 
cases,  the  interferometer  must  be  configured  such  that  orthogonal  SOPs  can  interfere.  There 
are  many  ways  to  achieve  this. 

A  first  way  is  to  put  an  analyzer  at  the  input  of  the  interferometer,  as  depicted  in  Figure  C.1. 
However,  if  no  polarizer  is  placed  at  input  and  both  arms  of  the  interferometer  have  no  effect 
on  the  SOPs,  no  cross-correlation  interferogram  representative  of  PMD  is  observed.  If  no 
polarizer  is  set  at  the  interferometer  input,  something  else  must  be  done. 

Second,  a  wave  plate  in  one  interferometer  arm  may  be  used  in  case  of  an  air-type  interfero- 
meter. Generally  speaking,  the  roundtrip  in  the  two  interferometer  arms  of  any  dual-path 
interferometer  can  be  represented  by  Jones  matrices  and  This  is  equivalent  to  a  wave 
plate  with  Jones  matrix  T  =  T^T2^  in  one  arm  only.  In  case  of  a  fibre  interferometer,  a 
Lefebvre  loop  may  be  put  in  one  arm,  and  adjusted  until  T  =  T^T2^  gives  the  desired  effect 
(a  given  cross-correlation-to-autocorrelation  ratio). 

One  particular  case  consists  of  putting  a  quarter  wave  plate  in  one  arm  of  a  Michelson 
interferometer  (or  a  half  wave  plate  in  one  arm  of  a  Mach-Zehnder  interferometer);  with  this 
configuration,  only  the  cross-correlation  interferogram  is  observed. 

C.1. 5      Polarization  scrambler 

In  Figure  C.2c,  the  polarization  scrambler  allows  the  selection  of  any  SOPs  for  the  FUT  input 
and  output.  The  polarization  beamsplitter  allows  simultaneous  detection  of  what  would  be 
detected  by  two  orthogonal  analyzer  settings.  The  functionality  of  the  polarization  scrambler 
selecting  various  SOPs  for  the  input  and  various  analyzer  settings  at  the  output  can  be 
achieved  by  other  means. 
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C.1.6      Polarization  beam  splitter 

A  polarization  beam  splitter  (PBS)  may  be  used  as  shown  in  Figure  C.2c  to  obtain 
interferograms  from  output  SOPs  that  are  orthogonal  (opposite  on  the  Poincare  sphere)  for 
the  same  l/O-SOP  combination.  These  two  interferograms  allow  the  calculation  of  the 
autocorrelation  and  cross-correlation  as  separate  functions.  Together  with  the  detection 
system,  the  PBS  forms  a  polarization  diversity  detection  system.  Means  other  than  the  PBS 
may  be  used  to  obtain  these  interferograms  from  orthogonal  output  SOPs. 

C.2  Procedure 
C.2.1  Calibration 

The  equipment  is  calibrated  by  checking  the  mechanics  of  the  delay  line  with  a  birefringent 
fibre  of  known  PMD  delay.  Alternatively,  an  assembly  of  birefringent  fibres  of  known 
characteristics  may  be  measured.  The  environment  and  the  launching  fibre  shall  be  stable 
during  the  measurement  period. 

C.2. 2      Routine  operation 

One  end  of  the  fibre  under  test  is  coupled  to  the  polarized  output  of  the  polarized  light  source. 
The  other  end  is  coupled  to  the  interferometer  input.  This  can  be  done  by  standard  fibre 
connectors,  splices  or  by  a  fibre  alignment  system.  If  the  latter  is  used,  some  index  matching 
oil  at  the  joints  avoids  reflections. 

The  optical  output  power  of  the  light  source  is  adjusted  to  a  reference  value  characteristic  for 
the  detection  system  used.  To  get  a  sufficient  fringe  contrast  the  optical  power  in  both  arms 
shall  be  almost  identical. 

C.2.2.1       Procedure  for  TINTY 

A  first  acquisition  is  made  by  moving  the  mirror  of  the  interferometer  arm  and  recording  the 
intensity  of  the  light.  The  fringe  pattern,  P{t)  ,  is  calculated  by  subtracting  the  "d.c."  part  from 

the  interferogram:  P{t)  =  P{t)- Pq  .  The  fringe  envelopes  that  are  generally  displayed  are  the 

absolute  value  of  the  fringe  pattern.  Typical  examples  of  fringe  envelopes  for  negligible  and 
random  polarization  mode-coupling  are  shown  in  Figure  0.3. 

In  case  of  insufficient  polarization  mode-coupling,  or  in  case  of  low  PMD,  it  is  recommended 
to  repeat  the  measurement  for  different  SOPs  or  to  modulate  the  SOP  during  the 
measurement  in  order  to  obtain  a  result  which  is  an  average  over  all  SOPs. 
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PMD  program  -  [PMD  Graph  :  c:\pmd\cablo102.pmd] 


'  I  Files     Parameters     Measurement     Graph     Options  ? 


PMD  Graph  :  c:>pmdlcablo102.pmd 
Delay  [pico-seconds] 
-5     -4     -3     -2      -10      1       2  3 
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PMD  Delai/  :  1.65  [ps] 
PMD  Coeff  :  1.126  [ps  /  Km-^1/2] 
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Figure  C.3a  -  Random  mode-coupling  using  a  TINTY-based  measurement 

system  with  one  I/O  SOP 


PMD  program  -  [PMD  Graph  :  c:\pmd\hibi2.pmd] 


■  I  Files     Parameters     Measurement    Graph     Options  ? 
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Figure  C.3b  -  Negligible  mode-coupling  using  a  TINTY-based 
measurement  system  with  one  I/O  SOP 

Figure  C.3  -  Fringe  envelopes  for  negligible  and  random  polarization  mode-coupling 

C.2.2.2       GINTY  procedure 

The  combination  of  a  particular  input  polarizer  setting  and  an  orthogonal  pair  of  analyzer 
settings  is  called  an  I/O  SOP.  Complete  the  scan(s)  for  the  two  interferograms,  from  the  two 

orthogonal  analyzer  settings  and  subtract  the  "d.c."  part  from  each  to  obtain  Px{t)  and  Py(r), 

the  orthogonally  generated  fringes. 
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The  cross-correlation  and  autocorrelation  fringe  envelopes,  EJt)  and  Eq(t),  are  calculated  as: 


(C.2) 


These  functions  are  squared  for  the  purposes  of  later  calculations  and  display.  Some  example 
squared  cross-correlation  results  are  shown  below.  Note  that  the  autocorrelation  peak  seen 
with  the  TINTY  is  not  present. 
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cr^  =  50  fs;  a  nearly-Gaussian  smoothed  envelope;  smoothing  is  for  guiding  the  eye  only:  analysis  is  not  performed 
on  any  kind  of  fit. 

Figure  C.4a  -  Random  mode-coupling  using  a  GINTY-based  measurement  system 

with  l/O-SOP  scrambling 
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Figure  C.4b  -  Negligible  mode-coupling  using  a  GINTY-based  measurement  system 

with  l/O-SOP  scrambling 
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'  -30      -24       -18       -12       -6         0  6  12        18         24  30 

Delay  ps 


NOTE  One  random-mode-coupling  section  with  Ul^.  =  10  +  one  negligible-mode-coupling  section  with 
DGD  =  P/WDRandom/4;  PMD  =  9,97  ps,  cr^  =  50  fs  with  nearly-flat  envelope. 

Figure  C.4c  -  Mixed  mode-coupling  using  a  GINTY-based  measurement  system 

with  l/O-SOP  scrambling 

Figure  C.4  -  Fringe  envelopes  for  negligible  and  random  polarization  mode-coupling 

(Ginty  procedure) 

Options  for  polarization  scrambling  follow.  When  multiple  I/O  SOPs  are  measured,  index  them 
with  /  for  later  calculations. 

C. 2. 2. 2.1        The  9-states  Mueller  set 

The  sum  of  nine  squared  envelopes  observed  with  nine  specific  I/O  SOPs  is  rigorously  equal 
to  the  uniformly  scrambled  mean  squared  envelope.  These  nine  I/O  SOPs  are:  three  analyser- 
axes  forming  a  right-angled  trihedron,  for  each  three  input  SOPs  also  forming  a  right-angled 
trihedron. 

C.2.2.2.2       Random  scrambling 

Scan-to-scan  scrambling:  automatic/manual  setting  of  the  scramblers  at  each  scan. 

Continuous  scrambling:  when  squared  envelopes  are  summed,  scrambling  can  be  performed 
while  scanning.  Automated  scramblers  are  set  to  cover  the  sphere  continuously  as  a  function 
of  time. 

Fast,  single-scan  scrambling:  if  scramblers  are  sufficiently  fast,  well-scrambled  squared 
envelopes  can  be  observed  in  a  single-scan.  However,  this  requires  special  provisions  to 
avoid  crosstalk  between  the  a.c.  part  and  the  previously  d.c.  part  of  the  interferogram. 

C.3  Calculations 

Both  these  calculations  report  PMD^^^.  For  ideal  random  mode  coupling  when  the  interfero- 
gram is  of  an  ideal  Gaussian  shape,  this  may  be  converted  to  PMO^^g  using  Equation  (3). 
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C.3.1      TINTY  calculations 

In  the  case  of  negligible  mode-coupling  the  PMD  delay  is  determined  in  the  same  way  for 
TINTY-based  and  GINTY-based  measurement  systems,  from  the  separation  of  the  two 
satellite  coherence  peaks  each  delayed  from  the  centre  by  the  DGD  of  the  fibre  under  test  as 
shown  in  Figure  C.Sb  and  Figure  C.4b.  For  this  case  the  DGD  is  equivalent  to  the  PMD  delay. 

<At>  =  —  (C.3) 

Co 

where  AL  is  the  moving  path  of  the  optical  delay  line  between  the  two  satellite  coherence 
peaks  and  Cgthe  light  velocity  in  free  space.  The  PMD  coefficient  is  given  by  <At>IL  where  L 
is  the  fibre  length  in  km. 

The  following  calculations  are  suitable  for  the  random  mode-coupling  regime  associated  with 
long  fibres/cables  or  links.  The  spread  in  the  fringe  envelope  -  discounting  the  centre  peak  - 
is  characterised. 

The  P/\//Dp|\/|3  value  is  determined  from  the  second  moment  (RMS  width)  of  the  cross- 
correlation  function  of  the  detected  signal  (fringe  envelope). 


PMDrms  =  J-  C7,  (C.4) 


where  Gs  is  the  RMS  width  of  the  cross-correlation  envelope. 

A  detailed  algorithm  for  the  calculation  of  cr^  from  a  measured  fringe  envelope  is  described  in 
Clause  D.1. 

For  certain  assumptions  given  below,  Equation  (C.4)  can  be  related  to  Equation  (3)  as: 

Ar2)  =  A(cTf)  (C.5) 

Equation  C.5  is  obtained  from  the  theory  given  the  following  assumptions: 

-  ideal  random  mode  coupling; 

NOTE  1  Ideal  random  coupling  means  l / h  ^ oo ,  and  a  uniformly  distributed  birefringence  axis.  L  is  the  device 
length  and  h  is  the  polarization  coupling-length.  For  a  device  consisting  of  N  concatenated  birefringent 
segments  of  length  h,  this  corresponds  toA/^oo  with  uniformly  distributed  axes. 

NOTE  2   Analysis  if  no  or  negligible  mode  coupling  is  possible. 

-  a  purely  Gaussian  source,  with  no  ripples; 

-  PMD  »  Oq,  where  Oq  is  the  RMS  width  of  the  auto-correlation  envelope; 

-  ergodic  conditions 

NOTE  3  Given  that  the  source  is  Gaussian,  the  result  must  be  some  form  of  weighted  average  of  the  DGD  values. 
This  weighting  is  not  specified  in  TINTY,  but  is  in  GINTY.  For  this  reason,  this  method  is  expected  to  give  a 
different  result  for  a  given  wavelength  range  and  time  than  one  of  the  methods  that  use  a  rectangular  weighting 
(such  as  for  example,  the  RTM).  The  assumption  of  ergodic  conditions  does  result  in  the  validity  of  the  expected 
value  relationship.  In  practice,  the  wavelength  ranges  sampled  by  different  implementations  of  the  other  methods 
will  also  vary,  which  would  imply  getting  different  results  between  them  as  well. 
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C.3.2      GINTY  calculations 

GINTY  allows  some  of  the  assumptions  required  for  Equation  (C.5)  with  TINTY  to  be  removed 
[4]: 

-  the  assumption  of  ideal  random  mode  coupling  is  not  required; 

-  the  assumption  of  a  Gaussian  source  is  not  required; 

-  the  assumption  that  the  PMD  is  large  compared  to  the  width  of  the  autocorrelation  function 
is  not  required. 

Form  the  mean  squared  cross-correlation  and  autocorrelation  envelopes,  Ex  (^")and  Eq{t)  as: 

/  / 

where  N  is  the  number  of  I/O  SOPs. 

Calculate  the  RMS  width  of  the  two  sampled  mean  squared  envelopes,  Oq  and  a^, 
respectively.  A  sample  algorithm  for  this  calculation  is  given  in  Clause  D.2.  The  mathematical 
definitions  of  these  widths  are: 


I  r2  (r)(E2  (r))  d  r  jr^  (r)(E2  (r))  d  r 


The  expected  value  operator  in  the  above  equations  is  with  respect  to  a  uniform  and  random 
sampling  of  the  l/O-SOPs. 

The  PMDp|\/|3  value  that  is  reported  is: 

n1/2 


PMDi 


RMS 


(C.8) 


Equation  (C.8)  relates  to  Equation  (3)  as,  from  [4]: 

(C.9) 


Using  the  ideal  definitions  of  the  RMS  width  terms  from  Equations  (C.7),  Equation  (C.9)  is 
exact  for  any  DGD  curve  present  at  the  time  of  measurement  and  any  source  spectral 
characteristics.  The  left  side  of  Equation  (C.9)  is  a  spectrally  weighted  (by  squared  power) 
RMS  calculation. 


Using  the  assumption  of  ergodic  conditions: 


'jAr2(v)s2(.)d, 
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Annex  D 

(informative) 

Determination  of  RIVIS  width  from  a  fringe  envelope 


This  annex  presents  two  methods  for  the  deternnination  of  the  RMS  width  fronn  fringe 
envelopes.  These  methods  are  mainly  used  for  the  interferometric  method  of  measuring  PMD. 
Clause  D.1  is  for  an  envelope  with  the  autocorrelation  peak  present  and  is  suitable  for  the 
TINTY  analysis.  Clause  D.2  is  for  envelopes  without  this  peak  and  is  suitable  for  the  GINTY 
analysis. 

D.1     RMS  calculation  for  TINTY 

Figure  D.1  shows  a  fringe  envelope  with  an  autocorrelation  peak  in  the  centre. 


G.650.2_FI.1 


Mirror  position  t  ps 


Figure  D.1  -  Parameters  for  interferogram  analysis 

Let  / J  denote  the  measured  intensity  of  the  fringe  envelope  at  increasing  positions  tj, 
j=  ^...N,  with  [tj]  =  ps. 

Step  1  -  Computation  of  the  zero  intensity  /o  and  the  noise  amplitude  A/a 
Definition:  A/5  =  round  (5  A//100) 


2N. 


(D.1) 


A/5  ^ 


;=1 


2N. 


(D.2) 


Na  =  ^X2-ll 
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(D.3) 


step  2  -  Definition  of  the  shifted  intensity  /, 


I J  :=  I J  -  Iq  if  I J  -Iq>  4Na 


/y:=0  iflj-lo<4Na 


Step  3  -  Computation  of  the  centre  C  of  the  interferogram 


IS/ 


c  = 


N 

M  

N 


step  4  -  Removal  of  the  central  autocorrelation  peak 

Definition:     y';/  =  the  largest  index  j  such  that  C  -  tj  >  (1-7) 
y'p  =  the  smallest  index  j  such  that  tj  -  C  >  (1-8) 

where     is  the  source  coherence  time. 

NOTE  1    For  cross-correlation  interferograms,  the  following  definition  shall  be  applied: 

J'r  :=y/+1 

Step  5  -  Computation  of  the  second  moment  S  of  the  interferogram 


S  =  - 
2 


Step  6  -  Truncate  the  interferogram 

Set  y^jn  to  the  largest  index  y  such  that     C  -  tj>  2S 

Set  y^ax  smallest  index  y  such  that  tj  -  C  >  2S 
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Step  7  -  Computation  of  the  second  moment  cr  of  the  truncated  interferogram 


y-ymin 


Jl 

y-ymin 


(D.11) 


{t-cf 

2 

Step  8  -  Computation  of  the  cr of  the  Gaussian  e  such  that 


Vmax 


{t-cf 


{t-Cfe    2cT^  df 


e  2fT^ 


(D.12) 


Step  9  -  Determination  of  P/WDr^s  <A^2>i/2 


2  ^1/2_  3 


(D.13) 


D.2     RMS  calculation  for  GINTY 


The  following  algorithm  results  in  a  robust  RMS  width  for  either  the  composite  square  cross- 
correlation  or  autocorrelation  envelopes  using  Method  C  (GINTY). 

The  algorithm  is  iterative.  For  a  given  iteration,  the  full  data  array  is  divided  into  two  sets:  the 
central  portion,  M,  containing  the  signal  and  the  tails,  T,  which  contain  noise.  Each  iteration 
results  in  different  definitions  for  these  sets.  The  result  converges  when  either  the  computed 
RMS  width  ceases  to  vary  or  when  the  set  definitions  stabilise.  For  an  iteration,  the  number  of 
data  points  in  each  set  are  denoted  as  nm  and  Nj. 

Let  Ij  denote  the  measured  intensity  of  the  envelope  at  increasing  positions  fy(ps),  (ps), 
J=  1...N. 

The  initial  definition  of  the  set  T  is  the  first  and  last  5  %  of  the  whole  array. 
Step  1:  Computation  of  the  zero  Iq 

l  =  ^'ljlNT  (D.14) 
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Step  2:  Definition  of  the  shifted  intensity  Ij 

lj=lj-lo  all  N  (D.15) 
Step  3:  Computation  of  the  centre  C  of  the  interferogram 

C  =  A^  (D.16) 

JeM 

step  4:  Computation  of  the  RMS  width  a  of  the  squared  envelope 

2  _  JeM  

JeM 

step  5:  Redefine  the  sets 

Define  M  as  the  set  of  points  for  which  C  -  4a  <  tj  <  C  +  4a. 
Define  T  as  the  rest  of  the  points. 

Step  6:  Repeat  steps  1  to  5  until  the  results  converge. 


(D.17) 
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Annex  E 

(informative) 

Glossary  of  symbols 


Symbology 

Arg 

ASE 

BBS 

C/Cq 

DGD 

DOP 
E 

EC 

f(AT) 

FA 
FT 

GINTY 
INTY 

I/O 

JME 

k 

Ic 
L 

LED 

N 

PbW 
PjojW 

P,{T)/P,{T) 

PBS 

PDL 

PDV 

PMD 

PMDq 

PSA 

R{A) 

RBW 

RTM 

s 

SOP 


Argument  function 

Annplified  spontaneous  ennission 

Broadband  source 

Velocity  of  light  in  vacuunn/in  free  space 

Differential  group  delay 

Maxinnunn  DGD  value 

Degree  of  polarization 

Number  of  extrema  in  R(A)  (Method  A) 

Extrema  counting 

Maxwell  probability  distribution 

Fixed  analyser  (Method  A) 

Fourier  transform 

General  analysis  for  Method  C 

Interferometry  method  (Method  C) 

Input/output 

Jones  matrix  eigenanalysis  (Method  B) 
Mode-coupling  factor 
Coupling  length 

Length  of  fibre/fibre  cable  test  sample 
Light  emitting  diode 

Total  number  of  measurements/population  of  mode-coupled 
fibres/wavelength  intervals 

Optical  power  recorded  with  analyser  in  place  (Method  A) 
Probability  of  exceeding  DGD^g^ 

Optical  power  recorded  with  analyser  rotated  90°  (Method  A) 

Optical  power  recorded  with  analyser  removed  (Method  A) 

Received  power  in  the  two  orthogonal  SOP  axes  corresponding  to  the 
fringes  in  Method  0 

Polarization  beam  splitter 

Polarization  dependent  loss 

Polarization  dispersion  vector 

Polarization  mode  dispersion 

Link  design  PMD  value 

Poincare  sphere  analysis  (Method  B) 

Output  ratio  from  PMD  measurement  system  (Method  A) 

Resolution  bandwidth 

Reference  test  method 

Normalized  output  Stokes  vectors 

State  of  polarization  (Method  B) 
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SPE  Stokes  parameter  evaluation  (Method  B) 

T  Jones  matrix 

Inverse  of  the  Jones  matrix 
fc  Optical  source  coherence  time  (Method  C) 

TINTY  Traditional  analysis  for  Method  C 

a  Single  parameter  which  specifies  a  Maxwell  distribution 

Chi-squared  variable 

Ah lAv lAq lAc    Finite  differences  computed  from  the  Stokes  vectors 

SA  Wavelength  step  size 

AA  Optical  source  spectral  width  (full  width-half  maximum  (FWHM)  unless 

otherwise  noted) 

Sv  Optical  frequency  step  size 

AO  Rotation  angle  on  Poincare  sphere 

St  Arrival  time  of  different  polarization  components  of  a  signal  or  pulse 
broadening 

^^%iax  Maximum  (Jrvalue  that  can  be  measured 

A^Jr^jp  Minimum  (^rvalue  that  can  be  measured 

At  DGD  value 

Ar^ax  Maximum  DGD 

<At>  Average  DGD  over  wavelength  scan  range  or  PMD  value 

<Az2>i/2  r.m.s.  DGD  over  wavelength  scan  range  or  PMD  value  (Method  0) 

<At>o  Maximum  PMD  specification  that  each  fibre  must  meet  in  a  population  of 

mode-coupled  fibres 

<At>f  Average  DGD  over  time 

<At>j  Average  DGD  over  temperature 

<At>ji  Average  DGD  over  wavelength 

Aco  Angular  frequency  variation  in  Method  B 

A  Test  wavelength  used  to  measure  PMD 

Aq  Central  wavelength  of  the  light  source 

A^IA2  First/last  wavelength  in  set  of  test  wavelengths  (or  position  of  first/last 
maximum  or  minimum  in  R(A)  in  Method  A) 

V  Optical  light  frequency 

p>\lp2  Complex  eigenvalues  of  T{a)+  Aco)T-\co) 

a  One-standard-deviation  uncertainty 

cTq  r.m.s.  width  of  the  squared  envelope  of  the  autocorrelation  interferogram 

(Method  C,  GINTY) 
ox  r.m.s.  width  of  the  autocorrelation  envelope  (Method  C) 

(Tp  Second  moment  of  FT  data  (Method  A) 

ax  r.m.s.  width  of  the  squared  envelope  of  the  cross-correlation  interferogram 

(Method  C,  GINTY) 

(jg  r.m.s.  width  of  the  cross-correlation  envelope  (Method  C,  TINTY) 

CD  Angular  optical  frequency 

Q  PDV 
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